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We use a simple elastic Hamiltonian for the vortex lattice in a weak impurity background which 
includes defects in the form of integer-valued fields to calculate the free energy of a vortex lattice 
in the deep H C 2 region. The phase diagram in this regime is obtained by applying the variational 
approach of Mezard and Parisi developed for random manifolds. We find a first-order line between 
the Bragg-glass and vortex-glass phase as a continuation of the melting line. In the liquid phase, 
we obtain an almost vertical third-order glass transition line near the critical temperature in the 
H — T plane. Furthermore, we find an almost vertical second-order phase transition line in the 
Bragg-glass as well as the vortex-glass phases which crosses the first-order Bragg-glass, vortex-glass 
transition line. We calculate the jump of the temperature derivate of the induction field across this 
second-order line as well as the entropy and magnetic field jumps across the first-order line. 



PACS numbers: 74.25.Qt, 74.72.Hs 

I. INTRODUCTION 

The phase diagram of high-T c superconductors as a 
function of the magnetic field H and temperature T is 
mainly governed by the interplay of thermal fluctuations 
and quenched disorder [l], Q , leading to various different 
states of the vortex matter summarized in Fig. 1. 

At low magnetic field the vortex solid melts into a 
vortex liquid (VL) via a first-order melting transition. 
Prominent examples of high-T c superconductors exhibit- 
ing a solid-liquid melting transition are the anisotropic 
compound YB^CuaOr-,? (YBCO), and the strongly lay- 
ered compound Bi 2 Sr 2 CaCu208 (BSCCO). The position 
of the melting line in the H — T plane is mainly influ- 
enced by the anisotropy of the superconductor In 
YBCO with its low anisotropy, most of the melting line 
B m lies in the vicinity of the upper critical field H C 2, 
i.e. B/H C 2 > 0.5 where B m H is the induction field 
in the relevant regime. The phase diagram for supercon- 
ductors with a melting line in this regime was discussed 
by us in Ref. 0|. This magnetic field regime consists of 
a vortex lattice, or vortex fluid, with overlapping vor- 
tex cores, if we use elastic constants calculated in Refs. 
H, @. The phase diagram was derived on the basis of a 
defect melting model set up in Refs. The model is 

Gaussian in the elastic strains and takes into account the 
defect degrees of freedom by integer- valued gauge fields. 
From this we derive effective Hamiltonians for the low- 
temperature solid and high-temperature liquid phases by 
summing over all defect fields. By further integrating out 
vortex degrees of freedom we obtain the partition func- 
tions of both phases. This is done with the help of the 
variational approach of Mezard and Parisi [i| , originally 
developed for random manifolds and applied later to vor- 
tex lattices without defects in Refs. [HJ, [U A similar ap- 
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FIG. 1: Sketch of the phase diagram of BSCCO or similar 
high-T c superconductors where the phase transition lines lie 
far below H C 2 . The solid line denotes a first-order phase tran- 
sition line being a first-order melting transition between the 
BG1-VL, BGl- VGl phase and a disorder induced first-order 
transition between the BG2-VG2 phase. The (blue) dashed 
curve denoted by Td is found to be a second-order glass tran- 
sition line The intersection point with the first-order line 
is denoted by GP2. The (red) dashed-dotted curve is the T x 
line found by Fuchs et al. using surface barrier experi- 
ments. It intersects the first-order line in the point GP1. 



proach was used recently to calculate the glass transition 
line for YBCO via Ginzburg-Landau theory [l2| . 

When including weak pinning, the solid phase becomes 
a quasi- long-range ordered Bragg-glass [1|. At higher 
magnetic fields, the quasi-long-range order is destroyed 
and there exist also a vortex-glass phase. The transition 
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is marked by the disappearance of Bragg peaks in scatter- 
ing data. We obtain in Ref. a phase diagram consisting 
of a unified first-order phase line between the Bragg-glass 
phase and the vortex-glass phase and the vortex-liquid 
which is sketched for BSCCO in Fig. 1. We point out 
that the first-order character of the transition line be- 
tween the Bragg-glass phase and the vortex-glass phase 
is not experimentally established yet for YBCO. It is de- 
duced from magnetic anomalies in response to the exter- 
nal magnetic field. For BSCCO, the first-order charac- 
ter of the Bragg-glass, vortex-glass line was measured by 
supercooling [13| and magnetic field shaking techniques 



Beside the unified first-order line found within our 
model for YBCO Q| (seen before within the Ginzburg- 
Landau model in Ref. fl5T) . a third-order glass transition 
line emanates near the critical point on the melting line 
as a phase boundary between the Bragg-glass and the 
vortex-liquid phase. We have shown in Ref. by us- 
ing hyperscaling relations that the higher-order charac- 
ter (more than second-order) of this line is in accordance 
with experiments and numerics which determines the 
scaling of the disorder phase correlation length [3, [TtJ • 
This glass transition line exists also for BSCCO. But be- 
side this parallel of the glass transition lines of YBCO 
and BSCCO, Beidenkopf et al. [3 found for BSCCO an 
additional second-order glass transition line in the Bragg- 
glass phase by using also the magnetic shaking technique. 
The line showed up by plotting the derivate of the mag- 
netic induction field with respect to the temperature. A 
jump was observed which also exist for the glass tran- 
sition line in the vortex-glass phases. Thus in contrast 
to YBCO, they found for BSCCO that both lines are of 
second-order characteristics. Both lines meet in a single 
point within experimental uncertainties. This point is 
not the critical point found for general doping [19 ] which 
is characterized by a vanishing of the entropy jump [l4[ 
being the maximum of the unified first-order line. We 
label both second-order lines in Figure 1 by Td- The in- 
tersection with the first-order line is denoted by GP2. 
Both lines divide the vortex-glass phase named by VG2 
from a phase named VG1 in Fig. 1 lying in the high- 
magnctic-field part above the first-order line. In the low- 
magnetic-field part the Td line divides two Bragg-glass 
phases denoted by BG2 and BG1 in Fig. 1 

Finally, we show in Fig. 1 a possible additional phase 
boundary labeled by T x which was found by Fuchs et al. 
[20| by measuring the vortex penetration through surface 
barriers. A similar line was also found by magnetization 
measurements [2lj . This line divides the vortex liquid 
denoted by VL from the phase VG1 shown in Fig. 1. 
It meets the melting line in a point to be referred as 
GP1. Note that it is not experimentally clear whether 
the T^-line has the characteristic of a phase transition. 
The X^-line does not correspond to the irreversibility line 
where magnetization sweeps show hysteresis. The posi- 
tion of this line in the case of BSCCO is mostly influ- 
enced by surface barriers 22j in contrast to YBCO where 



the pinning mechanism is responsible for the irreversibil- 
ity. This leads to a coincidence of the irreversibility line 
with the glass transition line between the vortex-glass 
and vortex-liquid [23|. It was shown in Ref. [2l[ via 
Josephson plasma experiments that the T^-line is not ac- 
companied by a Josephson decoupling between the layers 
ruling out the possibility of a transition from vortex lines 
to weakly coupled pancake vortices. That such a transi- 
tion exist was proposed theoretically in Refs. l24l. l25l, I2H 
So far we point out, that it is not experimentally clear 
yet what kind of phase VG1 is [2(|. There are, for exam- 
ple, hints that this phase could be a disordered flux line 
liquid [27[. This is suggested by muon spin rotation ex- 
periments which, however are in contradiction to other 
experiments which reported Bragg-peaks in this phase 
[28|]. In the interpretation of VG1 as a disordered flux 
line liquid, VG2 consists of a quasi-two dimensional vor- 
tex solid. 

It is the purpose of this paper to investigate the above 
phase transitions in the defect melting model mentioned 
earlier d, 0, [I| used in Ref. to calculate the phase di- 
agram of YBCO. We will first review briefly the model. 
A more comprehensive discussion can be found in our 
former papers and the book In contrast to YBCO, 
BSCCO is a strongly layered material where the coupling 
between the layers is described by the Josephson coupling 
in the Hamiltonian of the system. For high magnetic 
fields beyond the first-order line, one obtains a suppres- 
sion of the Josephson coupling between the layers [2l|, [29[ 
with respect to the electromagnetic coupling. In our elas- 
tic model with defects, we can not get this decoupling. 
We shall accommodate it effectively via an appropriate 
modification of the elastic moduli of the vortex lattice 
system in this region of the magnetic field. With the 
help of the elastic moduli of Brandt [|[ for BSCCO we 
show that one expects two thermal decouplings for the 
vortex strings in the liquid phase, corresponding to the 
two glass transition lines in Fig. 1, in contrast to the 
single line in YBCO [4]. On this way, we carry out the 
Mezard-Parisi analysis for the Hamiltonian of the vortex 
lattice system with pinning. It consists of a variational 
approach to fit the free energy of the replicated system 
with the free energy of a quadratic Hamiltonian. We ob- 
tain an almost vertical third-order depinning glass tran- 
sition line in the liquid high-temperature phase located 
in the vicinity of the T x line in Fig. 1 separating a full 
replica symmetric saddle point solution at high tempera- 
tures and a full replica symmetry broken solution at lower 
temperatures. We show that the saddle point equation 
to the variational free energy has no solution for very 
low temperatures. This is also the case when going be- 
yond the Mezard-Parisi theory within variational pertur- 
bation theory [30j | . This is a systematic extension of the 
Mezard-Parisi theory to higher orders. It is well-known 
phenomenon of higher-order variational perturbation ex- 
pansion of the quantum mechanical anharmonic oscilla- 
tor, that variations of the trial free energy do not neces- 
sarily have to show a minimum or a maximum [30| , where 
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the odd orders of variational perturbation theory possess 
a minimum, but even orders have no saddle points but 
only turning points. It is shown in Ref. [3l| for the anhar- 
monic oscillator that also turning points are acceptable. 
This is the principle of minimal sensitivity. Motivated by 
good results for the anharmonic oscillator we generalize 
the variational approach of Mezard-Parisi by using the 
principle of minimal sensitivity for the calculation of the 
variational free energy. With the help of this extension 
we shall obtain a variational free energy in the whole 
interesting regime of the H — T phase diagram for the 
vortex lattice. This phase diagram looks rather similar 
to the phase diagram in Fig. 1. The glass transition line 
Td corresponds then within our model to the tempera- 
ture where saddle point solutions to the variational free 
energy stops to exist and the best solution corresponds to 
turning points at lower temperatures where these turn- 
ing point solutions are still full replica symmetry broken. 
The transition show second-order characteristic and can 
be interpreted as a thermal depinning transition where an 
almost equally displaced substring due to disorder form- 
ing a plateau decouples from the impurities due to tem- 
perature fluctuations. 

Finally, we point out here as was also the assumption 
for YBCO [H that we will only consider the phase dia- 
gram in the H — T regime in the vicinity of the melting 
line. Going beyond this restriction would take much more 
effort being out of the scope of this work. 

The paper is organized as follows: In Section II and 
Section III we state the model and the Mezard-Parisi ap- 
proach to the free energy of the vortex lattice system for 
BSCCO. In Section IV we discuss solutions of the saddle 
point equation within the Mezard-Parisi approach. In 
Section V we consider the regimes where this equation is 
not solvable. Section VI goes beyond lowest order vari- 
ational perturbation theory using generalized principle 
of minimal sensitivity. Section VII discusses observable 
consequences of our theory. In the Appendices A and B 
we supply additional material to Sections V and VI. 



with the low-temperature Hamiltonian 

H [ui] = flr—oM = | y^(Vj^)(cn - 2c 66 )(ViUi) 

X 

+ 7}(^i u 3 + VjMi) C 66 (Vjltj + VjUi) + (V3W1) C44 (V3W1) 
= g ^2(^i u L) c ll (VjUl) + (V 3 U L ) C44 (V 3 U L ) 

X 

+ {ViU T ) c 66 (Vjitr) + (V 3 w T ) C44 (V 3 u T ) . (2) 

Here = P^u is the longitudinal part of the dis- 
placement where the projector Pl is given by (-Pi)jfe = 
-(l/i/|Vf|)Vj ® (l/vWiI)Vfc- The transversal part of 
the displacements is then given by uy = P^u = u u^. 
By using the dual representation in the form of stress 
fields we obtain a high-temperature representation of the 
partition function. This partition function describes the 
vortex system in the fluid phase. We obtain a partition 
function of the form {!]) with Hamiltonian 

H [ui\ = i? T ^ooK] = I X^ VlUl )( Cl1 _ c 66)(ViUi) 

X 

+ (V 3 Uj)c 44 (V 3 tij) 

v \ ~\ 

= ^ 2_^{ViU L ) (en - c 66 )(V 1 ul) + (V 3 m l ) C44 (V 3 u L ) 



+ (V 3 u T ) C44 (V 3 it r ) 



(3) 



and N = l/(4w/3) N . In the following, the subscripts 
i,j have values 1,2, and l,m,n have values 1,...,3 
where N denotes the number of lattice sites. The pa- 
rameter (3 is proportional to the inverse temperature, 
(3 = v c§(, / ksT (2tt) 2 , where the volume v of the funda- 
mental cell is equal to \fia 2 a^/2 for the triangular lat- 
tice. Here a is the transverse distance of neighboring 
vortex lines, and a 3 the persistence length of the disloca- 
tion lines introduced in Ref.[3|. Note that a 3 is assumed to 
be independent on the disorder potential in the average 
[32^ . Its value is given by [|[ 



II. MODEL 



The partition function to be used to describe the vor- 
tex lattice without disorder was proposed in Ref.0. It is 
motivated by similar melting models for two-dimensional 
square Q and triangular [1] crystals. In Ref. we derived 
from this partition function a low-temperature represen- 
tation. This corresponds to the partition function of the 
vortex-lattice in the crystalline phase 
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The lattice derivates Vi are built from the link differ- 
ences around a plaquette in the triangular lattice. These 
expressions can be found in Refs. 0, 0- By analogy V 3 is 
the lattice dcrivatc in z-dircction. 

The second term in the exponent of (fTj) 



H d i S [ui] = ^V(x + u), 



(5) 



accounts for disorder. We have suppressed the spatial ar- 
guments of the elastic parameters, which are functional 
matrices Cy (x,x') = Cy(x — x')._ Their precise forms 
were first calculated by Brandt [f| and generalized in 
Ref. [1] by taking into account thermal softening rele- 
vant for BSCCO. The elastic moduli C44 and cqq at low 
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reduced magnetic fields B/H C 2 < 0.25 are given by 
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l + A^g/(l + A^I z ) 
l + A 2 o6 fc 3 2 /(l + 2Ay(u2)) 



(7) 



where A c is the penetration depth in the xy-plane, and 
Kbz is the boundary of the circular Brillouin zone 
K 2 = AirB/fo. For BSCCO we use the two-fluid model 
H| A(T) = A(0)[1 - (T/T C )T 1/2 , t(T) = £(0)[1 - 
(T/T c ) 4 ] 1/2 /[l - (r/T c ) 2 ] . When calculating c 44 in © 
we have used a momentum cutoff in the two- vortex inter- 
action potential k < v / 2/(m 2 ) 1/ ' 2 , rather than the inverse 
of coherence length l/£ a fc as in Ref. 0. We shall not 
give here the explicit functional dependence of the elas- 
tic module c\\. This can be found in Ref. One can 
show that en S> C44 and Cg6 in the vicinity of the melt- 
ing line 3. This leads to the conclusion that one can 
neglect longitudinal fluctuations in comparison to trans- 
verse ones in the interesting regime Q • This will be done 
in the following. 

The last term in C44 ([7|) comes from the electromag- 
netic coupling between the layers. Its first two terms are 
due to the Josephson coupling (both terms arc vanishing 
for \ c /\ a b ^ 00) resulting in a vanishing of these terms 
in the case of zero Josephson coupling. It is possible to 
find approximations for C44 in ([7]) leading to tractable re- 
sults for the calculation of the free energy expressions of 
([1]) . In Ref. 13 we used the approximations 



B<t>o 



c i4 (k, k 3 ) 



32tt 2 A 2 (1+A 2 X 2 ) 

ab^ 1 ab BZ-' 

Bfa ln(l+2B\ 2 ab /<f> c 



32vr 2 A 4 ,fe 2 

ab o 



Cjl for h > ^ 

^ab 



(8) 

which arc justified for \k 3 \ < ~n / a 3 . Here, the Lindcmann 
parameter c 2 L = (w 2 )t->o/ '° 2 restricted to the transversal 
fluctuations is given by 
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where the average is taken with respect to the low- 
temperature Hamiltonian ([2]) without disorder represent- 
ing the elastic energy of the vortex lattice. The momen- 
tum integrations in |9]) run over the Brioullin zone of 
the vortex lattice whose volume is Vbz = (27r) 3 /u, as 
indicated by the subscript BZ. Kj is the Fourier trans- 
form of Vj @. Approximation ([5]) for C44 is correct in 
the regime i?7r 3 A 2 h /80o m (l/c|) < 1 which is valid in 
the vicinity of the melting line [3j. In this regime we 
obtain that c^k^k^) is dominated by the last term in 
§7} for |fc 3 | < 7r/a 3 . For higher magnetic fields than the 
disorder induced first order BG2-VG2 line (see Fig. 1) 
we have £?7r 3 A 2 h /8</>o ln(l/c 2 ) > 1 [|[ meaning that the 



first term in C44 ([7]) is dominated over the third term in 
the region k 3 ss ft/ 0,3- This implies that the approxima- 
tion ([5]) would result in a wrong approximation for the 
magnetic field regime above the BG2-VG2 line. We can 
see from ([2]) and ([3]) that the string tension C44 is not 
renormalized going from the vortex lattice to the vor- 
tex liquid. For deriving the full clastic constants 
one uses a quadratic approximation for the Josephson 
coupling cosine phase difference term in the Ginzburg- 
Landau model for BSCCO. It was shown in Refs. l34l. l35l 
theoretically and in Refs. [2]], [2^ by determing Josephson 
plasma frequencies for BSCCO that one gets a suppres- 
sion of the full Josephson energy between the layers when 
going from the vortex solid to the vortex liquid crossing 
the BG2-VG2, BG1-VG1 line. This leads effectively to a 
softening of the Josephson terms of C44 in the VG1 and 
VG2 phases being the first two summands in ([7]). This 
justifies to use ([5]) as a good approximation for the full 
string tension in the whole interesting regime when also 
including Josephson decoupling. 

The disorder potential V^(x) due to pinning is assumed 
to possess the Gaussian short-scale correlation function 



y(x)y(x') 



d(T) a 3 



(10) 
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x'i) S : 



where K(xi—x' i ) l/(£') 2 for |x— x'| < and zero else- 
where. The parameter </>o is the magnetic flux quantum 
4>o = hc/2e, and parameter £' is the correlation length of 
the impurity potential, which has a similar value as the 
coherence length £ ab in the xy-plane. In the following, 
we use an effective disorder correlation function with the 
Fourier transform 



A"(<7) = 2™xp(-£'\ 2 /2) 



(11) 



leading also to an exponentially vanishing of the disorder 
correlation function in real space. In Ref.i we have used 
this form for the correlation function in the solid phase 
for YBCO. In the present material BSCCO, this is even 
more justified because the disorder potential looks <5-like 
for the vortices, due to the large thermal fluctuations of 
the vortices near the melting transition line 0]. 

The temperature dependence of the parameter d(T) 
has two sources. One is the temperature dependence of 
the correlation length, the other is based on the pinning 
mechanism where we discuss in the following the ST C - 
pinning or Si -pinning mechanisms [l[ 

d(T) = doil-T/Tc)- 1 / 2 for ST C - pinning ,(12) 
d{T) = d {l-T/T c ) 3/2 for pinning. (13) 



III. MEZARD-PARISI METHOD 

We now carry out the calculation of the partition func- 
tion ([T]) which is still complicated due to disorder. In 
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Ref. we have done this for YBCO by using a quadratic 
approximation in the disorder strength. This leads to 
a reentrant behaviour of the melting line in the H — T 
plane which did not agree with experimental results. By 
using the variational approach of Mezard-Parisi Q to go 
beyond the quadratic approximation this reentrant be- 
haviour is disappeared, leading to good results for the 
form of the melting line and agreement to the transition 
line between the Bragg-glass and vortex-glass. Here we 
use again the Mezard-Parisi theory to perform a sim- 
ilar calculation in the case of BSCCO. In order to go 
beyond second-order perturbation theory in the impu- 
rity potential, we use first the well known replica trick 
[361 ]. The Mezard-Parisi theory consists in replacing 
the non-quadratic part of this replicated Hamiltonian as 
quadratic with a possible mixing of replica fields. By us- 
ing the Bogoliubov variational principle we can find the 
best quadratic Hamiltonian so that its free energy named 
F var is as close as possible to the actual free energy of the 
system. This means that we have to search the minimum 
of 



Fyax — Atrial + (H — -fftrial)trial 

with the harmonic trial Hamiltonian 
v 
2 



fftr 



ial 



£5>«(x)G-i(x-xV(x'). 

X,X' CKj/3 



(14) 
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Here (-)trial stands for the average with respect of the 
Gibb's measure of the trial Hamiltonian iftriai, while H 
denotes the replicated Hamiltonian. The indices a, (3 
denotes the replicas. 

In the general form, the search for an extremum is 
a complicated problem. A strong simplification for this 
was founded by Parisi for random-spin systems where 
he suggest to deal with a trial Hamiltonian within some 
subalgcbra known as the Parisi algebra. This restriction 
can be motivated by physical arguments (37[. It will be 
clear soon for the solid as well as the fluid phase that the 
transverse part of G a p can be chosen to have the form 



G a l — G Q 5 a f} + (Tap 



(16) 



where Go is the transverse part of the Green function of 
the Hamiltonian Ho[ui] ([2]) in the solid phase and Q in 
the fluid phase. 

Within the Parisi-algebra, the self-energy matrix a a p 
depends effectively only on one parameter [§] (see also 
App. B). In the general form it is allowed to be a con- 
tinuous function cr(s) with < s < 1 Then the 
variational free energy has the form 0, 0] 



A/ var = AF vaI /N = / V ar(B[A]) - / var (0) 
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N is the number of lattice sites, and G 1 is given by 
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in the solid low-temperature phase corresponding to lp|). 
and 



Go '(Ms) 
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(21) 



in the liquid high-temperature phase corresponding to 
Q. Here e; are the three unit link vectors around a 
plaquette in the triangular lattice. The gap function A(s) 
and the self-energy function cr(s) are related by 



A(.>-/ ^ 
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B[A(s)] is given by 
B[A(s)] = ^-g[A(s)] - M! f ds>±g[A(s>)] 

V S V J s s z 

k B T 



knT 



ff[A(l)] 



ds'a'(s')g'[A(s')} . (23) 



In order to find a saddle point of F val /N we have to 
take the derivative of (fT7| with respect to A(s). This 
results in ujl 



a(s) 



k B T 



V (2B[A(s)]) 



(24) 



where T>'(x) is the derivative (d/ dx)T>(x). The disorder 
function T> is given by Q 
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(k B Tf X\ b e + (u 2 ) 



(25) 



In the following, we discuss solutions of this equation 
in the cases that <r(s) does not break the replica symme- 
try, possesses one-step replica symmetry breaking, or a 
continuous replica symmetry breaking. 

In order to solve (|24|) . we first have to calculate g(A) 
(fT9| which we will denote by g T ^°(A) with ([20]) in the 
solid phase, and by g T ^°°(A) with JSJ) in the fluid 
phase. We shall use the elastic constants cqq of Eq. (J6]) 
and the approximation (|8|) for C44. In the liquid case, the 
result is 
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and for the vortex solid 

T~>0fA\ ~ 0.0987T aa 3 V3 1 a 2 ^1/2 



9 



J (A) 



+1.38- 



1 



,(2) 
44 



(l + Z z (0) A/2) 



(27) 



where A = Aa\/c^, and C44 denotes the function 
C44(k, fc 3 — * 0) of Eq. ([8]) for & 3 < l/\ a b, and c 44 denotes 
C44 of Eq. ([8|) in the region fc 3 > 1/X a b for ^3 = 1/03, i.e. 
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For the derivation of (|27p we have used the approximation 

o^cee/a 2 cfJ -C 1 valid in the vicinity of the melting line. 
We used further the abbreviation 
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(29) 



In the solid phase the following abbreviations will be use- 
ful 



HO) 



Z 



(i) 




»1, 



(30) 



1.38^— 



bib*) »i(31) 
a a 



Note that A a f, « a in the vicinity of the critical point on 
the melting line where the disorder is most infhicntial the 
shape of this line. In this regime we obtain large numbers 
on the right hand sides of ([29]). (|30j). and (01]). In (26} 
and (|27|) the last terms have their origin in the integration 
(fl"9f over momenta 1/A a b < |&3| < 7r/a 3 . The other terms 
come from the integration over small momenta. Expres- 
sions g T ^°(A) d27j) and g T ^°°(A) ([26]) are not exact 



results of the integration in Q19p . They are good approx- 
imations for g(A) but also for <?'(A) and <?"(A) in the re- 
gion A < 1 in the fluid phase and A < ((A a6 /a) 2 /Z ; (0) ) 2 / 3 
for the solid. It will be seen below that these are the rel- 
evant regimes for F var . 

We now define the quantity 



A = 



k B T 



<l:>, 



(32) 



which will be useful below. Comparing g (A) in ([2"B"]) and 
(|27|) with the corresponding expressions for YBCO we 
obtain that only the last terms are different. The first 
term in ([2"6"| leads in the case of YBCO to the decou- 
pling scenario in the fluid high-temperature phase [l[. 



This is the regime where temperature fluctuations starts 
to dominate over disorder fluctuations for the coherently 
pinned vortex line pieces given by T>(0)A ~ 1. The length 
of such line pieces are given the Larkin length L c where 
disorder fluctuations grow to value 

One can now show by generalizing the calculation of 
the vortex fluctuations due to pinning and thermal fluc- 
tuations for YBCO in Rcf. [38| that the additional last 
term in ([2"6"]l causes a new length scale beyond the Larkin 
length. At this length scale, the vortex fluctuations are 
approximately constant forming a plateau. This length 



scales like Lt 



a^zf*^) for thermal fluctuations, 

for disorder fluctuations at low 
temperatures. Beyond these lengths, the displacement 
fluctuations starts to increase proportional to the cubic 
distance due to pinning, and proportional to the distance 
for thermal fluctuations. Both lengths, that of the co- 
herently pinned vortex line pieces but also the vortex 
substring on the plateau can decouple due to thermal 
fluctuations. 

Below we find two different depinning phase transition 
temperatures: One takes place when the temperature 
fluctuations exceed the disorder fluctuations over the co- 
herently pinned vortex line pieces where the the Larkin 
length fulfills L c > Lt- This leads to the well-known 
depinning temperature of the coherently pinned vortex 
substring given by T>(0)A ~ 1 corresponding to the third- 
order phase transition (|56|) below. The second depinning 
transition takes place when the temperature fluctuations 
exceed the disorder fluctuations over the plateau in the 
regime where the Larkin length is given by L c < Ljj. 
This leads to the depinning temperature T>(Q)A ~ zf^ 1 
corresponding to the second-order phase transition (|57|) 
given below. Of course one can also see both depinning 
temperatures mentioned above in the temperature de- 
pendency of the Larkin length L C (T) (39|, 



and 



like L D = a 3 (^ (0) ) 2/3 



IV. SOLUTION OF THE MEZARD-PARISI 
SADDLE POINT EQUATIONS 

In the following, wc discuss the solutions of the 
Mezard-Parisi Eqs. (J22J), (J23J) and ([24|) in the liquid and 
the solid phase. 



A. Liquid Phase 

In order to solve the Mezard-Parisi equations we trans- 
fer the analysis for YBCO of Ref.ito BSCCO. Note that 
by neglecting the second term in ([27)) we obtain a similar 
expression for 5 T ^°°(A) as for YBCO 0. This leads 
to the following results: The stable solution for A(,s) 
is replica symmetric for V(0)A < 2/V% and full replica 
symmetry broken in the case V(0)A > 2/V3- That for 
example the one-step replica symmetry breaking solution 
is not stable can be seen from the following fact: The 
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one-step self-energy function Ai is given by 



V(0) 



(33) 



Ai 



log 1 



(l + ^4>) 



AO) Ai 



1+Z 



(0) Aj 



V3 ^4 



Aj 1+Z 



(0) Aj 



= 1 



The one-step symmetry breaking solution of the saddle 
point equation (|24p is stable when the replicon eigenvalue 



A = 1- 



8 



Ai +4 



2?(0) 



(34) 



^(0) A X 



(i +z <°»^; 



i + -2-^- 



7 (<>] 



'Ai l + Z\ 



(0) A x 



is larger than zero. By comparing (|33|) with ([34]) we 
obtain that A < when Z, (0) Ai/2 < 13 meaning that the 
one-step replica symmetry breaking solution is unstable 
in this range. More generally one can show similarly as 
in Ref. that all finite step replica symmetry breaking 
solutions are unstable for V(0)A > 2/^/3. 

Thus, we expect a continuous replica symmetry break- 
ing solution in this parameter range. Note that continu- 



ous step replica symmetry breaking solutions of the sad- 
dle point equation ([24]) are stable in general [3, S3| ■ We 
now calculate this solution by using the methods given 
in Ref. 0. First, the full replica symmetric solution for 
V(0)A < 2/V3 is given by 



AT 

A(s) = for ^-T>(0)A < 1 . (35) 



The continuous replica symmetry broken solution for 
T>(0)A > 2/V3 can be derived from the saddle point 
equation ([24]) by differentiating both sides with respect 
to s resulting in 



ct'(s) = -a 1 { 8 )4 



g'[A(s)]V" (2B[A(s)}) . (36) 



This means that tr(s) is either constant or solves equation 
(|36jl divided by cr'(s). Dividing (|36p by cr'(s) and forming 
the derivate with respect to s, we obtain with 



g'[A{ S )] 2 V" (2B[A{s)\) 



-sg"[A{s)]V"{2B[A{ s )}) . 



(37) 



Equations ([3"B"]) and ([57]) can be solved algebraically for 
the unknown functions A(s) and -B[A(s)] leading to 



A(s) = 



5(0) Xl/2> 



1+Z\ ' A 



1+(Z 1 (0, ) 2 A 3 / 2 



for 

/ r- \V 3 

f ^X>(0)A) s for 
(^D(O)A) for 



s < 



(^27(0)A)V3 



— 7= i < S < S c , 

(^•D(0)A)V3 - 



< s < 1 



(38) 



where we used Z\ /A ss c 2 L a 2 /2£ 2 ^ 1 in the vicinity of 
the melting line Q and the abbreviation zf^ = zf' /(1+ 
A/2). This means that A is constant for s > s c . The 
constant s c is given by the equation 

(i + (zr) 2 A3/ 2(Sc) ) 5/3 fv - 3 y/3 

^ — — '- =\l±V^)A\ s c . (39) 

l + |(^°))3A^( Sc ) V 2 / 

Finally, we can calculate the disorder part of the vari- 
ational free energy A/ var (fT7]) in the liquid phase [!]• 
With A/ var = A/J£ + Af$£ where A/*» is the kinetic 
part represented by the first term in (|17[) of the disor- 
der energy and A/J^j.* is the potential energy part of the 



disorder energy (second term in (|17p ). we obtain 
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and 



2/3 



glass transition line separating the phases VG1-VL at 
T>(0)A = y/3/2 is of third-order. We found the same 
order for the depinning transition in YBCO 0]. 



ds- 



A 1 /2( s ) 



(^■O(O)A) 1 / 3 



_(l + (Z i ( ° ) )2A3/2( S )) 

1/3 \ / 



1/3 



+ (1-Sc) 



A 1/2 (s c ) 



l + (Z ; (0) ) 2 A 3 /2( Sc ) 



1/3 



(41) 



By taking into account (f38|) and (f39|) we obtain that the 



B. Solid Phase 

In the solid phase one can show that finite-step replica 
symmetry breaking solutions arc unstable Similar 
as in the discussion of the fluid phase in the last subsec- 
tion wc obtain the following continuous replica symmetry 
broken solution of the saddle point equation (|24| 



A(.s) 



A 



l+2Z i (0, ^ 1) A3/ 2 3( 

A1/2 .098 3 (1+ZW) 3 
1 + Z^A 1 / 2 



C6sa 3 
64tF I r ( 1 )„2 

I — 



1/2 



2/3 



S for < S < S c , 

(I?(0)A) for s c < s < 1 



(42) 



where we used zf" 1 /A ps c|a 2 /2£' 2 > 1 in the vicinity of and 
the melting line [3 and the abbreviation =Z^/(1 + 
zf^ K/2). The constant s c is given by equation 



A 2 / 3 ( Sc ) 



(i + z<pK^{ Sc f 



l + 2Z ; (0) zi 1) A 1 /2( Sc ) 
(2^) 2 / 3 



3(V3/2) 1 /3 



(P(O)A) 



1/3 ( c 66«3 



2/3 



(43) 




1-s 



(i+zi i} Ai/2( S )y 

AV6( Sc ) 



1/3 



(V(O)A) 2 / 3 



(45) 



V. EXISTENCE AND STABILITY OF SADDLE 
POINT SOLUTIONS 



Finally, wc can calculate the disorder part of the vari- 
ational free energy A/ var p7|) . We obtain for this energy 
in the solid phase 



A f kin = - 

J var o 1 r I 2 

2 l07T \ C66fl3 



x/3 k B T 



Jds± ( 2 -A^ 2 ( s) + ^A 2 ( s ) 




A 3 / 2 (.s c ) 



,(i) 



-A 2 ( Sc ) 



(44) 



Trying to solve the implicit equation for A(s c ) in the 
liquid phase (last line in (|38| ) and the solid phase |42|) 
we obtain that in both cases a solution is not existing 
for very large V(0)A > (2?(0)A) max corresponding to 
low temperatures or large disorder strengths according 
to (J25J), (f32]l. 

We obtain from Eq. (|38|) taken at s = s c or directly 
from (|36|) for (2?(0)A) max at low temperatures 



(O(0)A) n 
A B 

(S c )n 



Z 7 (0) 

3 2 / 5 

(^(°))6/B 



A, 



(46) 

(47) 
(48) 
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where A = A max and s c = (s c ) max at T>(0)A 
(D(0)A) max . In the solid phase we have from (j42jl 



(V(0)A)„ 
A,. 

(Sc)n 



,(0) 



9 

5 ' 



(49) 
(50) 

(51) 



The calculation was done by maximizing the left hand 
side of the implicit equation of ([55)1 and (|4"2"j) for s c < 
s < 1 with respect to A. This then gives the maximal 
2?(0)A value given by (T>(Q)A) max where wc still get a 
solution for both implicit equations. Summarizing we 
obtain that the continuous replica symmetry broken so- 
lutions in the liquid as well as the solid phase stops to 
exist for (X>(0)A) max w 2Z/ 0) /V3. It was shown in the 
last paragraph that stable solutions of the saddle point 
equations are infinite replica symmetry broken (where 
we have the restriction to V(0)A > 2/\/3 in the liquid 
phase). More generally we obtain that every saddle point 
solution of (|24p irrespective of its form is unstable for 



V3T>(0)A/2 > Z\ because the replicon eigenvalue 



A = 1 



g'(A(l))V" 



2^fl(A(l)) 



(52) 



where A(l) is the self-energy function at s = 1, is nega- 
tive in this range. 

We point out that in the liquid phase similar to the 
continuous replica symmetry breaking solution discussed 
above also the one-step replica symmetry breaking solu- 
tion ([33]) is no longer existing for V3V(0)A/2 > 4Z, (0) . 
This can be seen by Taylor expanding the left hand side 



of 



33]) with respect to Z^Ai/2. Note the difference 
in the prefactor of zf ] compared to (|4"6j). ((49|) . This 
leads us to the more generally assumption that there is 
no saddle point solution of (f2~4"|) for \fZT> (0)^4/2 larger 

than rsj Z'i ■ . This is proved in Appendix A where it is 
shown that this is true for every finite-step replica sym- 
metry broken solution of (|27|) by using results derived in 
Ref. |. 

We point out that T>(0)A ~ Z^ is in fact a relevant 
parameter region for the glass transition line because we 
expect that the critical point is around T>(0)A ~ z\ a \ 
Here we use (|2T|) with the fact that the quadratic ap- 
proximation to the disorder energy at the peak should be 
approximately k B T 0], i.e. V{c 2 L a 2 )k B T ~ k B T. Note, 
as is shown in Fig. 1, the glass transition lines separating 
the phases BG2-BG1, VG2-VG1 crosses the first-order 
line BG2-VG2. BG1-VG1 in the vicinity of the critical 
point for optimal doping. 

In the variational perturbation treatment of the anhar- 
monic quantum mechanical oscillator we obtain a similar 
phenomenon. The even variational approximations to 



the free energy posses no extremum in the variational 
parameter [30( . Only the odd perturbative orders where 
the Mezad-Parisi theory belongs to the lowest order ap- 
proximation within this perturbation theory, has a true 
minimum. In order to see whether we have a similar 
situation here, i.e. whether higher order variational ap- 
proximations to the free energy posses a physical plausi- 
ble extremum for A(s), we will calculate in the following 
section higher order variational approximations to the 
free energy. 



VI. BEYOND LOWEST ORDER VARIATIONAL 
PERTURBATION THEORY 

In this section, we will go beyond lowest order varia- 
tional perturbation theory outlined in Section III. Start- 
ing from (|14p we can immediately write down the next 
beyond lowest orders of the free energy F within varia- 
tional perturbation theory [30l ] 



F v: 



ial 



(A - A tr 



ial/ trial, c 



m w 

En 



1=2 



11 (k B Tf 



— {{H — Atrial) )trial,. 



(53) 



((A - #triai) J )triai,c is the the averaging of (A - Atrial)' 
with respect to the trial Hamiltonian (|15p where we only 
take the connected expectation value part which means 
for example in second-order ((A — trial) 2 ) trial, c = ((A — 
iftrial) 2 )trial - ((A - fftriai)) 2 ria i- In order to calculate 
the free energy i^ var as in Section III within m'th order 
variational perturbation theory we limit the Z-sum in (|53p 
to I = to. F val . corresponds to the exact free energy of 
the system for to — > oo which means that F V&T does not 
depend on the choice of the trial Hamiltonian Atrial- The 
truncated sum depends on the choice of Atrial- Since the 
infinite sum is Atrial independent, the best truncated sum 
should depend minimally on A tr i a i. A first approximation 
would be in taking a saddle point of F V3-r with respect to 
the trial Hamiltonian At r i a i leading to ([24]) in the case 
to = 1. 

To calculate F var beyond lowest order for a trial Hamil- 
tonian At r i a i within the Parisi algebra is not an easy 
task. When going beyond lowest order we expect that 
the continuous replica symmetry breaking self-energy 
functions are still the relevant ones as was shown in 
Section IV via stability considerations for the Mezard- 
Parisi case corresponding to first order variational per- 
turbation theory. We carry out the calculation of the 
free energy in Appendix B within second-order varia- 
tional perturbation theory (m = 2). We show that for 



,(o) 



there exist also in this case 



V(0)A > (V(0)A) max 
no continuous solutions of the saddle point equation in 
this second-order case. Thus, in contrast to the anhar- 
monic oscillator where the variational perturbation the- 
ory leads to a solution of the saddle point equations for 
every odd order (30l ] a similar phenomenon is not exis- 
tent in our case. Without explicit proof we now state the 
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conjecture that this is true for every finite order within 
variational perturbation theory. This means that there 
exist no saddle point of F var for large V3V(0)A/2 > zf 1 ' 
which is a relevant physical regime outlined at the end 
of the last section. One way out of this dilemma is to 
continue the continuous replica symmetry broken solu- 
tions given in (|38p for the liquid and (|4"2"]l for the solid 
to the regime V(0)A > (V{0)A) max w 2Zf ] /y/3 by look- 
ing closer to the anharmonic oscillator problem solved 
in Ref. [50 via variational perturbation theory. As men- 
tioned above, for equal orders within the variational per- 
turbation expansion which means m € 2Z in (|53p one 
does not find a saddle point with respect to the trial har- 
monic Hamiltonian -fftriai being a quadratic potential in 
the anharmonic oscillator case. There it is shown that 
one gets good accordance with numerical solutions of the 
Schrodinger equation when interpreting the requirement 
of the minimally dependence of Fva, r on .fftriai mentioned 
above by a vanishing of the second-order derivation of 
F var with respect to fftriai- This is equivalent to the de- 
mand that the first order variation of -F V ar on -fftriai is 
minimal. 

Transforming this general rule to our case by using 
()24p . the self-energy function <r(s) is given by the mini- 
mum of the functional 



Min 



1 + 



2 k B T 
a(s) v 



V'(2B[A(s)]) 



(54) 



where we assume as a first approximation that this mini- 
mum is not dependent on s. This leads us to the following 
result for T>(0)A > (X>(0)A) max : 

The solutions A(s) of ([54]) where a(s) and A(s) related 
by are given by ([55]). |@2J) with the substitution 
T>(0)A — ► (D(0)A) max . The variational energies A/ var 
are given by (|40"|) . (|4Tj) in the liquid phase and (|44l . 
(|45| in the solid phase with the same substitution. Fur- 
thermore, one has to multiply the potential part of the 
disorder energies (|4"Tj) and (|45|) by a correction factor 
{V(0)A)/(V(0)A) max for V(0)A > (£>(0)A) max . 

Summarizing, we obtain for BSCCO a third-order glass 
transition in the liquid phase having its reason in the 
breaking of the full replica symmetry across the transi- 
tion line at T>(0)A = 2/V3. A similar transition was 
also found for YBCO Beside this transition we will 
show in the next section that (f5"4")) leads additionally to a 
second-order glass transition line at V(0)A w 2Zf^ j\fi 
in both phases. We point out that this transition is not 
reasoned in the generalization of the saddle point cri- 
tcrium for the variational free energy to the more gen- 
eral principle of minimal sensitivity (|54[) . Up to now, 
we have only searched a saddle point of the variational 
free energy in the self-energy matrices a a p of the Parisi 
form (see the discussion below (Tl"5"|) ') which could be mo- 
tivated physically [13] as an Ansatz for the glassy-state 
self-energy matrices. Nevertheless, it could also be pos- 
sible that the restriction to this subspace is the reason 
that we do not find a saddle point of the variational free 



energy for V(0)A > 2Zf ) jyf%. On the other hand it is 
clear that also in this case the leaving of the stable saddle 
point solutions from the subspace of self-energy matrices 
of the Parisi form leads in general to a non-analytically 

of the free energy at the point D(0)A w 2Zf^ jyf% and 
thus to a phase transition. 

As we explained above the reason that the saddle point 
solutions of the variational free energy stops to exist 
within the Mezard-Paris theory lies in the non-solvability 
of (|5rj|) for s = s c . This follows further from the fact 
that 5(A) ~ 1/A for large A and that V" (2B [A(s c )]) w 
X>"(0)£' 6 /ff [A(s c )] 3 for the relevant A(s c ) values where g 
begin to show the behaviour 5(A(s c )) ~ 1/A(s c ). In de- 
riving the approximation for T>" above we use c|et 2 3> £' 2 
for BSCCO (see the notes below (|3SJ) and i[g2"])t In con- 
trast to this we find for YBCO c|a 2 <C £' 2 leading to the 
existence of the saddle point solutions of the variational 
free energy in the whole H — T plane although we have 
also 5(A) ~ 1/A for large A in this case [4j. This is 
the reason that one does not find the second-order glass 
transition line in YBCO. 

Finally, we note that Giamarchi et al. in Ref. [l(| 
only consider the small A behaviour of 5(A) which is 
presumably the reason that they did not find the second- 
order glass transition line at least in the solid phase. The 
reason that we can compare only our low-temperature 
solid phase results with results in this paper lies in the 
fact that they did not consider defects as we do here being 
relevant in the high-temperature liquid phase. Note that 
they did not use temperature softened elastic constants 
in their calculation relevant for BSCCO [|[. 



VII. OBSERVABLE CONSEQUENCES 

In the following, we use the intersection criterium Q 
with variational free energies (TlTf and (fl8|) to get the 
first-order line separating the phases BG2-VG2, BG1- 
VG1, BG1-VL. This results in 



B m (T) 



1 



(1 - (T/T c ) 4 ) 2 0; 



x cxp 



192 v^tt 7 
2 



\ 2 



k R T 



(A/, 



T->-0 
var 



(0)A 2 (0) 
A/, 



{k B Tf 



T^oo\ 
var J 



(55) 



where Af^*°° is given by the disorder part of the vari- 
ational free energy which is the sum of (|40| and (|41 [) in 
the liquid case. A/^7^ corresponds to the disorder part 
of the variational free energy in the solid case given by 
gH and 05]). 

Beside this first-order transition line we obtained a 
third-order glass transition line of the depinning form 
in the fluid phase separating the VL and VG1 phases 



(56) 



and a second-order glass transition line separating BG1 
with BG2 in the solid phase and VG1 with VG2 in the 



liquid phase given by (|46| . (|49)) 

V(0)A = (P(0)A) raax « ±=Z\ 0) . (57) 

This means that we obtain within our analytical approx- 
imation a unified glass transition line in both phases in 
correspondence to the experimental findings shown in 
Figure 1. 

In the following figures, we use parameter values for op- 
timal doped BSCCO given by A q6 (0) w 2300A, £ a6 (0) « 
30A, Cu02 double layer spacing a s = 14A, T c = 90K and 
the anisotropy parameter 7 = A c / A a b ~ 250. Due to the 
small coupling between the layers, the Josephson form of 
the interlayer coupling leads to a non-negligible softening 
of A c or 7 = Xc/Xab, respectively, as a function of B and 
T. In Ref. [2^ it was found by Josephson plasma experi- 
ments that X C (B,T) is nearly of the form l/\ 2 c (B,T) ps 
{l + F{B I B m )) I T) with some function F which can 
be found in (2<| and further that Aj?(0, T)/\ 2 c (B m , T) tn 
0.6. This leads to 7 = \ c (B,T)/\ ab w 250 in the vicin- 
ity of the first-order line separating the phases BG2-VG2, 
BG1-VG1 BG1-VL, Here we used A c (0, T)/\ ab (T) « 200 
as in Ref. y. 

In Figure 2 we show (|55|) corresponding to the first 
order line separating phases BG2-VG2, BG1-VG1, BG1- 
VL for ST C pinning given by the correlation function 
(|12p (upper figure) and 81 pinning (|13[) (lower figure) for 
various constants do- The square points in the figure 
denotes the experimentally determined first-order BG2- 
VG2J3G1-VG1, BG1-VL line of Beidenkopf et al. in 
Ref. [l8|. The do values of the the straight (black) curves 
are chosen in such a way, that we reproduce in one of 
the best ways the experimentally given curves of Bei- 
denkopf et al. and also the glass intersection point GP2. 
We obtain do = 2.5 ■ 10~ 6 in the ST C pinning case and 
do = 4 • 10~ 6 for Si pinning. The curves of representative 
variations of these almost best parameter values are given 
by the (red) dotted lines in Fig. 2. We obtain discrepan- 
cies in the form of the first-order BG2-VG2, BG1-VG1, 
BG1-VL line from the experiment. There are a large vari- 
ety of the concrete forms of this line in the literature (see 
for example [l4| for an almost horizontal BG2-VG2 line 
with a small kink near the intersection point GP2). The 
reason for the discrepancies comes mainly from the sen- 
sitivity of the curve on the disorder function [1] but also 
the neglection of the layerdness of BSCCO in our case 
could be one factor. Without taking into account disloca- 
tions, we expect a Josephson decoupling transition which 
is nearly temperature independent for low temperatures 
[35l ] . The melting line and the decoupling line lies on 
top of each other when taking into account dislocations 
leading to the first-order BG2-VG2, BG1-VG1, BG1-VL 
transition line. Note that the Josephson decoupling is not 
complete over the transition line and further that the la- 
tent heat due to the Josephson degree of freedom is only 
16 % of the total latent heat over the first-order transi- 
tion line p9j . The competition between the temperature 
independent decoupling transition and the temperature 
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FIG. 2: The BG2-VG2, BG1-VG1, BG1-VL first-order tran- 
sition lines B m (T) given in (|55[1 for 5T c -pinning (upper fig- 
ure) and <5/-pinning (lower figure). The solid (black) lines are 
calculated with parameters for do which gives one of the best 
fits to the experimentally determined phase diagram in Ref.flSl 
(square points) within the pinning mechanism (do = 2.5 TO -6 
for <5T c -pinning, do = 4 • 10™ 6 for 51 -pinning). Dotted (red) 
curves are calculated by a variation of these parameters given 
by d^ = (1 ± 1/2) do- The vertical markers denote the in- 
tersection points of the VG2-VG1, BG2-BG1 glass transition 
line and the disorder induced first-order line BG2-VG2, BG1- 
VG1 named GP2. GPl denotes the intersection point of the 
third-order VG1-VL glass transition line with the first-order 
melting line BG1-VG1, BG1-VL. 



dependent three dimensional first-order line should take 
into account the correct form of the whole first-order line 
for layered materials. 

The small vertical marks on the curves in Fig. 2 de- 
notes the glass intersection point GP2. We obtain espe- 
cially for the ST C pinning case differences in the location 
of the glass intersection point GP2 with the experiment. 
In all shown three ST C pinning cases the glass transition 
point GP2 lies in the vicinity of the critical temperature 
T c where in both pinning mechanisms also the glass in- 
tersection point GPl is located. Summarizing, we obtain 
as was also the case for YBCO Q that the SI pinning 
mechanism gives a better accordance to the experimental 
curves and glass intersection points than the ST C pinning 
mechanism. 

In Figure 3 we show for do = 4 • 10~ 6 with SI pinning 
correlated impurities the whole phase diagram calculated 
with (j5"5"]) , (fB"6"|) and ([57)) corresponding to the parameter 
values do of the (black) solid line in the lower picture 




FIG. 3: Phase diagram for BSCCO in lowest order variational 
perturbation theory corresponding to the Mezard-Parisi the- 
ory. Lines represent the theoretical determined phase transi- 
tions between the various phases calculated for SI pinning with 
do = 4 ■ 10~ 6 corresponding the solid line in the lower picture 
in Fig. 2. Points represent the experimentally determined 
phase diagram of Beidenkopf et al. pjjfl. The solid (black) 
line denotes the first-order BG2-VG2, BG1-VG1, BG1-VL 
line calculated by (|55p . The glass transition lines BG2-BG1, 
VG2-VG1 are given by the (blue) dashed lines calculated from 
the numerical generalization to the approximation (|57p as de- 
scribed in Section V. The red dashed-dotted line is the glass 
transition line VG1-VL derived by using expression (f56)l 




FIG. 4: Entropy jump ASi per double layer and vortex (up- 
per figure), and the magnetic induction jumps AB (lower fig- 
ure) over the first-order line BG2-VG2, BG1-VG1, BG1-VL. 
Calculations are based on expressions derived in Ref. |j. In 
both figures we used for the (black) solid lines the ^-pinning 
mechanism with do = 4-10 -6 corresponding to the parameters 
of the phase diagram in Fig 3. The points in both figures rep- 
resent experimental values determined by Zeldov et al. [411 ]. 



in Figure 2. Note that H C 2 can not be resolved in this 
figure being almost vertical directed on the right bound- 
ary. Again we show for comparison the experimentally 
determined phase diagram of Beidenkopf et al. where 
the square points denote the first-order BG2-VG2, BG1- 
VG1, BG1-VL transition line. The (blue) circle points 
denote the experimentally determined T<j line BG2-BG1, 
VG2-VG1 of Beidenkopf et al.. This line has to be com- 
pared with the (blue) dashed lines VG2-VG1, BG2-BG1, 
calculated with |57|) where we get small discrepancies in 
the intersection point on the first-order line of the up- 
per high magnetic curve VG2-VG1 in the liquid phase 
with the small magnetic field transition line BG2-BG1 
in the solid phase. The reason is that we did not use 
the analytical approximation \/3(V (0)A) max /2 f=a zf' 
for (2?(0)A) max valid in both phases but the numerical 
determined values calculated from the condition that (|36|) 
stops to be solvable as described in Section V. As men- 
tioned by Beidenkopf et al. in Ref. [IH it could be exper- 
imentally possible that both lines do not intersect. From 
Fig. 3 we obtain that the point GP2 does not coinci- 
dence with the maximum of the theoretical determined 
first-order BG2-VG2, BG1-VG1, BG1-VL transition line 
which coincides with the critical point Q (see also Fig. 4). 
This is possible for general doping [l9j . Nevertheless, we 
obtain a discrepancy between the position of our glass 
transition lines and the experimental findings. One rea- 



son comes from the approximations to the elastic moduli 
carried out in Section II but also corrections to ([26| , |27|) 

where we used a 3 CQQ/a 2 c^ 1. These approximations 
getting worse for higher magnetic fields Q. This leads 
to an additional bending of the first-order line in the di- 
rection to the temperature axis shown in Fig. 1 of Ref. [H 
without pinning. To get the same effective bending of 
this line as in the experiments we have to use a smaller 
d value leading to BG2-BG1, VG2-VG1 lines located 
at smaller temperatures according to (|5T|) . Furthermore, 
a source of the additional bending can be also due to 
the decoupling transition between the Josephson layers 
as discussed above. 

Beside these reasons also the restriction to the low- 
est order variational perturbation approximation could 
be a source for the difference of our theoretical finding of 
the glass transition line and the experimental ones. The 
calculation of the free energy within second-order vari- 
ational perturbation theory is outlined in Appendix B. 
We did not carry out the calculation of the phase dia- 
gram within this order which is rather non-trivial being 
out of the scope of this work. 

Finally, the (red) dashed-dotted line in Figure 3 shows 
the VG1-VL glass transition line calculated by the help 
of the depinning temperature formula (|56p . We do not 



show for comparison the T x line of Fuchs et al. [20| in the 
figure because they did not use an optimal doped crystal 
in the experiment. 

In Fig. 4 we show in the upper picture the entropy 
jumps per double layer and vortex ASi and in the lower 
picture the magnetic induction jumps AB over the first- 
order BG2-VG2, BG1-VG1, BG1-VL transition line. The 
(black) full line is calculated with do = 4 • 10~ 6 in the 51 
pinning case corresponding to parameter values of the 
phase diagram in Fig. 3. We used formulas derived in 
Ref. for the calculation. Note further, as was also the 
case in Ref. Q, that we did not use corrections for AB 
by considering explicitly the difference of the induction 
field B and the external magnetic field H. These differ- 
ences are negligible in the interesting regime [Ilj]. The 
square points (blue) are experimentally determined val- 
ues measured by Zeldov et al. [4l[ for optimal doped 
BSCCO crystals. We note that there are other exper- 
iments in the literature for non-optimal doped crystals 
where AS and AB varies significantly [3, E2|- The rea- 
son for this difference is not clear. The largest differ- 
ence in Figure 4 between experimentally and theoreti- 
cally determined curves is at high-temperatures near T c . 
As noted in 0] this comes mainly from contributions of 
thermally activated vortex loops not inherent in our vor- 
tex lattice picture. 

In the paper of Bcidenkopf et al. [ijj the order of 
the glass transition lines VG2-VG1, BG2-BG1 was de- 
termined by measuring the magnetic induction field and 
its derivate with respect to the temperature across this 
line. They found a jump of dB(H, T)/dT across the line 
leading to the conclusion that this transition is of second- 
order. They also deduced from their experiment that the 
jumps over the glass transition line are of almost the 
same magnitude in the BGTBG2 phase and the VGT 
VG2 phase [HI, [l^] . Nevertheless the displayed curves in 
their paper show a much smoother behaviour of the mag- 
netic induction curve and its temperature derivate near 
the glass transition line in the BG2-BG1 solid phase than 
in the VG2-VG1 phase. The problem of determing the 
order of the transition comes mainly from a large noise 
on the magnetic induction curves having its reason pre- 
sumably in the spatial and temporal inhomogeneities of 
the s yst em. This is the reason that Bcidenkopf et al. in 
Ref. [l8| did not get a clear jump in the derivative in all 
measurements (while the induction itself bends sharply) 
to allow a systematic quantitative study of it 43| . 

The magnetic induction field B is given by 



B 



H 



B 



(58) 



where Z n n is the partition function of the static non- 
fluctuating part of the vortex lattice. In Fig. 5 we show 
the disorder part of the magnetic induction field (solid 





T[K] 



FIG. 5: The disorder part Sdia — otT of the magnetic induction 
field defined by (|59|l (solid (black) curves), and its derivate 
with respect to T, i.e. dBdis/dT — a (dashed (blue) curves) for 
two different magnetic fields B either in the VG2-VG1 phase 
(upper figure) or in the BG2-BG1 phase (lower figure). Here a 
is some subtraction parameter determined such that the £?dis 
curve is symmetric around the glass intersection temperature. 
We used a = 0.0055Gauss/A' in the VG2-VG1 phase and 
a = 0.009Gauss/i^ in the BG2-BG1 phase. The left-hand 
y-axis denotes the scale for _Bdis — otT, the right-hand y-axis 
for dB dis /dT - a. 



(black) curves) B given by 



d 1 

5 dis = -47r— -A/ v 
oB v 



(59) 



and also its derivate with respect to T (dashed blue 
curves), i.e. dBdis/dT as a function of tempera- 
ture for two different magnetic fields. As mentioned 
above, the magnetic induction contribution from the 
non-fluctuating part of the partition function log(Z n fl) 
is negligible in comparison to H in ([58]) to the induc- 
tion field B. This means that we obtain for the jump 
values over the glass transition line AdB(H,T)/dT 
AdB dis (B,T)/dT. 

The upper picture in Fig. 5 shows the disorder part 
of the magnetic induction field B^is in the liquid phase 
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for B = 550 Gauss. In the lower picture we show the 
disorder part of the magnetic induction field -Bdis in the 
solid phase for B = 180 Gauss. In correspondence to 
Beidenkopf et al. [l8T ] we subtract to B^is a term linear 
in the temperature T to get a symmetrical curve around 
the glass transition temperature. We obtain a negligi- 
ble jump AdBdis{B, T)/dT over the glass transition line 
in the solid BG2-BG1 phase. This is in contrast to the 
jump AdB<a s (B,T)/dT in the liquid high-temperature 
VG2-VG1 phase. By comparing the absolute values of 
this jump with the corresponding experimentally deter- 
mined jump values AdB(H, T)/dT determined in Ref.Qjl 
our values are about one order of magnitude smaller. 
Note that for our theory AdB dis (B, T)/dT is about 10" 1 
smaller in the BG2-BG1 phase than in the VG2-VG1 
phase. This could not resolved within our numerics in 
Fig. 5. That this is true can be seen from the following 
scaling consideration 



8 47T 

A-B dis (B,T) = --A 







5 d d 

"dBlAjs) A/var ar {sH dfJK(s 



A/ var A A ( s ) 



dsds' 
4tt 



TBi 



-(a/ v ^ + a/^M^ (0) a( Sc )) : 



(60) 



where we used ([54]) in order to substitute terms contain- 
ing A/p^* to terms containing A/™?. One can then see 
from analytic approximations but also numerical consid- 
erations that both terms in (|60p are of almost equal value 
in the liquid phase VG2-VG1 but that the first kinetic 
term of the disorder free energy in (|60p is much larger in 
the solid BG2-BG1 phase than the potential second part. 
Our numerics gives that the kinetic part of the disorder 
free energy A/ var in the BG2-BG1 phase is one order of 
magnitude smaller than in the VG2-VG1 phase. 

One source of the difference between the jump values 
of our theory and the experimental numbers could beside 
the approximations wc used in our theory also the addi- 
tional in-plane ac equilibrizing magnetic shaking field in 
the experiment of Beidenkopf et al. [43|. This shaking 
field is of the same magnitude as the dc magnetic field 
in z-direction. It is immediately clear from the results in 
Ref. [H[ as well as the theoretically and experimentally 
determined results for an additional in-plane dc field in- 
stead of the ac field H, Ef| that the shaking field 
has only a small effect on the position of the first-order 
line and also the jump values ASi and AB. This can 
be understand by using the anisotropic scaling theory [l[ 
leading to an attenuation of the in-plane field by a fac- 
tor A a h/A c . In contrast to this we obtain from Fig. 5 
that due to the smallncss of the magnetic field -Bdis the 
shaking field can still have an effect on the jump value 
AdB/dT across the glass transition line BG2-BG1, VG2- 
VG1. Note that an in- plane magnetic dc field can even 



put additional dislocations in the vortex lattice [46J, |4j 



VIII. SUMMARY 

In this paper, we have derived the phase diagram for 
superconductors which have their phase transition lines 
at magnetic fields much smaller than _ff c2 , i.e. B/H c2 < 
0.25 such as BSSCO. The model consists of the clastic 
degrees of freedom of the vortices with additional defect 
fields describing the defect degrees of freedom of the vor- 
tex lattice in the most simple way. For the impurity po- 
tential we have restricted ourselves to weak pinning 5T C 
and (5Z-correlated impurities [l[ . This model was formerly 
used by us for describing the phase diagram of supercon- 
ductors with a melting line near H C 2 4|. The layered 
structure of the superconductor, i.e. the Josephson cou- 
pling form between the layers, is not explicitly consid- 
ered. We take this special coupling only into account via 
the elastic moduli of the lattice and an experimen tally 
and analytically based decoupling scenario [U, [25|, [261 ] . 
In order to treat the impurity potential approximately 
we use a theory developed first by Mezard and Parisi 
Q for random-manifolds. This is based on a variational 
approach to the free energy via a quadratic trial Hamil- 
tonian. After stating our model in Section II we have 
discussed the Mezard-Parisi theory of the vortex lattice 
system in Section III. The minimum requirement for the 
trial free energy of the quadratic Hamiltonian leads to 
the saddle point equation (|24|) where the stable solu- 
tions are full replica symmetric for V(0)A < s/3/2 in the 
fluid phase with T>{0)A is defined in {55} and {32]). Ev- 
erywhere else, the solutions arc continuous replica sym- 
metry broken. We expand these solutions to low tem- 
peratures T*(0)A > (X>(0)A) raax with {ST}. That the 
saddle point equation (f24|) has no solution in general 
for V(0)A > (V(0)A) max is shown in Appendix A. The 
Mezard-Parisi theory is the lowest-order approximation 
of a more general perturbation theory known as varia- 
tional perturbation theory. In Appendix B we show how 
to go beyond the lowest-order approximation for the vor- 
tex lattice system up to second-order where also in this 
case a saddle point solution is not existent. Motivated by 
good results of the variational perturbation treatment for 
the anharmonic oscillator we generalize in Section VI the 
minimum requirement of the variational free energy de- 
terming the trial Hamiltonian to a more generalized prin- 
ciple of minimal sensitivity given in {54]) . This leads to a 
second-order phase transition line located at the points 
in the H — T plane where the saddle point solutions cease 
to exist. 

As was discussed by us at the end of Section VI, for 
YBCO in contrast to BSCCO the magnitude of the lat- 
tice fluctuations near the melting line is smaller than 
the correlation length of the impurity potential, i.e. 
(cLa ) 2 <C C 2 H- This is the reason that one does not 
find a similar non-existence of saddle point solutions to 
the variational free energy, in certain regions of the H — T 
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plane for YBCO as we obtain for BSCCO. This leads to 
the absence of the second-order phase transition line in 
YBCO. Further we note that Giamarchi and Doussal fl"ot | . 
who calculated the physics of the vortex lattice with pin- 
ning but without defects valid in the solid phase of real 
systems, did not find in their work the ceasing of saddle 
point solutions to the variational free energy in certain 
regions in the solid phase. The reason lies in the fact that 
they only consider small trial dimcnsionlcss gap functions 
A(s) in their calculation. Also they did not use temper- 
ature softened elastic constants relevant for BSCCO Q . 

The procedure described above leads to the following 
physical consequences for BSCCO. Due to the form of the 
clastic moduli in the deep H C 2 region, we obtain two glass 
phase transitions of the depinning form. The first transi- 
tion line of third-order is located in the fluid phase at high 
temperatures not far from T c . It is given by (|56jl identi- 
fied as the depinning temperature of a coherently pinned 
vortex substring. It separates the full replica symmetric 
solution to the variational energy at high temperatures 
(VL phase) and a full replica symmetry broken solution 
at lower temperatures (VG1 phase). This transition cor- 
responds to the glass transition in YBCO. The transition 
line is located in the vicinity of the experimentally found 
T x line [2(| . The second transition is of second-order ([ST]) 
dividing the Bragg-glass and the vortex-glass phases in 
four regions. It separates a full replica symmetry broken 
saddle point solution of the variational free energy (VG1, 
BG1 phases) and a full replica symmetry broken turning 
point solution (VG2, BG2 phases). This transition line 
is a temperature depinning transition where a substring 
which is almost equally displaced due to disorder forming 
a plateau decouples from the impurities due to temper- 
ature fluctuations. We find that the derivate jump of 
the magnetic induction field with respect to the temper- 
ature over this glass transition line in the Bragg-glass 
phase is negligible in comparison to the jump in the liq- 
uid phase. We compare this line with the experimentally 
found second-order glass transition line by Beidenkopf et 
al. [l8| located in the vicinity of our line. The jumps 
of the temperature derivate of the magnetic induction 
field in the vortex-glass phase of our theory is about one 
order of magnitude lower in comparison to the experi- 
mental values of Beidenkopf et al. [l|[. They obtain a 
similar value for the jumps in the Bragg-glass phase over 
the glass transition line. In comparison to the glass tran- 
sition line separating the VG2-VG1 phases they found a 
softening of the jump in the vicinity of the glass transition 
line in the BG2-BG1 phase consistent with our findings. 

We calculated the first-order melting transition line 
and its disorder induced continuation dividing the Bragg- 
glass phase BG2 and the vortex-glass phase VG2 by us- 
ing an intersection criterium for the low and the high- 
temperature expansion of the free energy. The whole 
theoretical determined phase diagram and the experi- 
mentally ones determined by Beidenkopf et al. [18[ is 
shown in Fig. 3. Finally we compared the entropy jumps 
per layer and vortex, and also the magnetic induction 



jump over the first-order line with the experimental find- 
ings of Zeldov al. [4l|. This is shown in Fig. 4. 

Summarizing, we have calculated the phase diagram of 
a vortex lattice model stated in 0, Q for BSCCO with- 
out taking explicitly into account the layered structure 
of the material. Although we found certain quantita- 
tive differences in the position of the experimental deter- 
mined phase transition lines, the overall phase diagrams 
looks rather similar. Discrepancies are maybe due to the 
approximative evaluation of the theory and the layered 
structure of BSCCO. 

We would like to thank H. Beidenkopf, E. H. Brandt, 
and A. Sudb0 for useful discussions. This work was sup- 
ported by Deutsche Forschungsgcmeinschaft under grant 
KL 256/42-2. 



APPENDIX A: GENERAL PROOF OF THE 
NON-EXISTENT OF FINITE-STEP SADDLE 
POINT SOLUTIONS FOR LOW 
TEMPERATURES 

In this section we show that there exist no finite-step 
saddle point solution for the variational free energy A/ var 
in the range T>(0)A > zf^ within first-order variational 
perturbation theory . This was shown in the continuous 
and additionally in the one-step case for the liquid phase 
in Section V. In order to derive this we use results derived 
in Section C in Ref. 0. We obtain for an i?-step replica 
symmetry breaking solution 



(Ef =1 



V(2^g[A mR }) 

where we used that Ao = 0, rriR+i = 1 and 

f A(s) d 

It is shown in Ref. that Z > 0. We have 
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in the liquid case and 



5(A) 



V3 1 (c^aA |2 A3/2 
2 8tt 1 c 66 a| I \ 3 



4Z, 



(i) 



log l + ^=| 



r(P) A 
'I ~2 



A 



(A4) 



16 



for the solid. Next we use the inequalities (X^i a i) 2 — 
i? X^ili a 2 for real number Qj, . . . ,or and further that 
A,. < A i+1 , Af/A 2 R < S*/Sx, (1/rm - l/m l+1 f < 



1/mj — l/mf +1 resulting in 



S(A mR ) 
RAl lR 



< 



V(2^fg[A 



■ 2 • ( A5 ) 

This inequality can be only fulfilled for y/3D(0)A/2 < 
ARZ\ 0) which is a generalization of the one-step replica 
symmetry breaking case in the fluid phase discussed be- 
low Eq. ([521). 



APPENDIX B: SECOND-ORDER VARIATIONAL 
PERTURBATION EXPANSION 

In this Section we calculate the second-order expansion 
terms within variational perturbation theory (|53|) . The 
aim is to show that also to this order there are no sad- 
dle points of Fvar ([53]) with D(0)A > zj ' corresponding 



to (|4"6")) and (|4"5)) in the first-order case. We restrict our- 
selves to solutions of the saddle point equation with full 
replica symmetry which were relevant in the first-order 
case according to Section IV. 

In the following, we will calculate first the expectation 
value of the disorder part of the replica Hamiltonian Q 



a,/3 



(Bl) 

in which a, (3 run over the replica indices. We allow only 
for onsite interactions which were justified in Ref. for 
YBCO. For the present compound this approximation is 
even more appropriate for BSCCO since £ a b ~ £' <C a. 
With this disorder part (|B1[) we obtain for H — i?triai m 
(El 

H - ff tria i = H dis - ~J2 E «r(x)o- a/ jt4(x) . (B2) 
x Q ,/3 

We now classify terms of higher-order variational pertur- 
bation theory in two groups. When expanding ((H — 
-Atrial)') trial hi (|53p we obtain first terms of the pure dis- 
order Hamilton form ((i?dis) z )triai,c which we denote by 
((H— J?triai) }triai,c l- Second, there are monomials which 
contain at least one self-energy matrix factor <j a p in it de- 
noted by {{H — i?triai)')triai,c,2- These terms can be most 
easily treated by the square root trick [3(| ■ We now cal- 
culate first terms of the pure disorder Hamilton form. 
1. Pure disorder terms in Hamiltonian 



As before, we restrict ourselves to the transver- 
sal part of the 2x2 Green function G(x) defined 
by 1/(2tt) 3 / d 2 qdq 3 (q T <g> qr) G{q)e i ^ +ic ' 3X \ 
The replica sum in (|B3[) is of the form 
J2 a 0y8 F[G a p, GjSt G ai , Gps, G a s, G ( 3 7 ] where F[-\ 
is some functional of the various Green functions. 
Since G a p is some matrix within the Parisi alg ebra 
the functional F has the ultrametric property [37l. l48j. 
Following Temesvari et al. j4^|, we denote the size 
of the Parisi blocks with p r , r = 1 . . . R, where R is 
the maximum level of replica symmetry breaking. We 
fix po = n and pr+\ = 1, the latter being the size of 
diagonal elements and n is the number of replica fields. 
The matrix elements a a j3 , that belong to the rth level of 
replica symmetry breaking arc all equal to some number 
<7 r , r = 0, . . . , R. The replica overlap function is defined 
by a n P — r when a a p = a r . 



The fact that the Green function G a p is in the Parisi 
algebra implies that the Green function definitely de- 
pends on the overlap a fl /3 which we denote in the fol- 
lowing by Ganp. Furthermore, the operation a fl /3 on 
the replica indices has the ultrametric property. This 
means that whenever we may choose three replicas a, 
/?, 7, cither all three of their overlaps are the same, i.e. 
an/3 = a:n7 = /3n7, or one e.g. a n (3 is larger than 
the other two. In the latter case the two are equal, i.e., 
af)P>aC)j — (3Dj. This means that of the three 
Green functions Gap, G a7 , and G/3 7 only two are dif- 
ferent. Similarly, of the six Green functions in F only 
three arc different. The various possible Green function 
combinations can be most easily determined by mapping 
these six Green functions onto the edges of a tetrahedron 
where the Green functions on the adjacent edges of a face 
must fulfill the ultrametric property. 



Within second-order variational perturbation theory, 
we obtain 

trial 
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In the following, we restrict us to the leading term 
x = x' in (|B3|) . By carrying out the q, q' integral we 
obtain 



x e 



■V(G 77 (0)+G 44 (0)-G 7i (0)-Gj 7 (0)) 
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((#dis) 2 )trial,c ~ (k B T) 2 ^ J2 P2 (°) 



(^f idaa ~ 9 a p) + (977 " 57*) + ^) J ' 
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The last subtracted term in (|B4j) is due to the connect- being the same equation as in the first-order variational 
edness of ((ifdis) 2 )trial,c- From (|B4|) we obtain perturbation theory case We now restrict (|B4[) to 

the Parisi algebra by carrying out the program outlined 
V" a a(i = , (B5) above leading for n -> to 
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(B6) 



Here g s is the momentum integrated Green function of 
G s according to Eq. (p~9|) . We define <5(s;— ~) by the 
functional J dsj<5(sj— = if (<?) where -ff is some 

functional of the integrated Green function g Si and g = 

For calculating the saddle point equation up to second 
order variational perturbation theory corresponding to 
(f24|) the derivate (i5/<5<7 s )((i?dis) 2 )triai,c is relevant which 
should be added with an appropriate factor to the right 
hand side of Eq. ([24]) . In order to derive this equation 
we first give the variational free energy F vav /N within 
second-order variational perturbation theory denoted by 

/var,2 



fvar.2 



2{k B T) 



{ trial j 



(B7) 



where / va r.i corresponds to the variational energy within 
first-order variational perturbation theory given in (|17p 
and (HHJ), i.e. / var> i = / var (0) + A/ var; i where A/ Vari i is 
a modification of A/ var specified in (fTT|) according to 



Pi A/, 



kin 

var 



pot 
var 



(B8) 



The additional prefactors Pj are modifications due to 
second-order perturbational expansion from terms pro- 
portional tO ((H — -fftrial) 2 



/trial,! 



B2|) containing at least 
one factor a a p. The constants Pi lying between 1/2 and 
3/2 will be determined in the next subsection. 

Carrying out the variation of / va x,2 with respect to g s 
we obtain 



P 3 cr(s) 



P 4 V'(2B[A( S )})- 



6 ((H dis ) 2 ) 

trial,* 

Sg~ s {k B Tf 



(B9) 

The calculation of ((5/(5g s )((_ffdis) 2 )triai,c is tedious but 
straight-forward. Due to lack of space, we do not state 
the result here. 

In order to discuss the sign of ((-ffdis) 2 )triai,c and 
(<5/<5p s )((-f/dis) 2 )triai,c we repeat the form of the Green 
functions in the Parisi algebra 0] 



1 



g - g s = 



(27T) 



ds'- 



d 2 kdk 3 
a'(s) 



G - 1 (k,fc 3 ) + A(l) 



(G - 1 (k,fc 3 ) + A( S ))^ 



(BIO) 
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and 



(2tt) 3 



1 



d 2 kdk 3 G (k, ha) 



ds- 



1 



A(s) 



,o S 2 G - 1 (k,fc 3 ) + A( s ) 
By using (|B7j) . jB9|) we obtain 



(Bll) 



s = s c for small temperatures. More precisely we find 
that ([Bi~2|) is not solvable for V(0)A > zf ] by using 
(V"(2B[A( Sc )])/v)- 1 D Sc (S/Sg Sc )((H dis f) triahc /(k B T) 3 ~ 

r»(o)A/z ; (0) . 



2. Terms containing at least one factor o 



a/3 



P 3 a'( S ) = -a'( S ) g'(A(s)) 



(B12) 



x ("^p 4 ^ ( 2 B[ A( s )]) +J D s / «^ is ) 2 )t 1 ,ai,c 



corresponding to (|36|) in the first-order case. Here, we 
have used the differential operator 

D s = (c/(s) g'[A(s)]\ '1. (B13) 



Dividing (|B12[) by <x'(s) and forming the derivate with 
respect to s, we obtain 



2k B T\ 2 g'[A( s )} 3 (2k B T 



g"[A{s)] v v 

-,2 $ ((-ffdis) 2 )trial, 



PJ)'" (2B[A(s)]) 



{k B Tf 



(B14) 



In contrast to the first-order results (|36|) , (|37|) the second- 
order variational perturbation equations (|B12p . (|B14|) 
give no longer local algebraic equations for i?[A(s)] and 
A(s) but integral equations involving both quantities for 
different s. 

From (|B6[) we obtain that i5/<5g<i((-ffdis) 2 )triai.c depends 
through g s on s. One can show after a tedious but 
straight-forward analysis that 



trial, c > o, 

6 



((tfdisrWc > o, 



-Ds 3 T ((-ffdis) 2 )trial,c > 0. 



(B15) 

(B16) 
(B17) 



Here s c is defined such that cr'(s) = for s c < s < 1. As 
in the first-order case, s c can be determined by (|B12|) 
for ,s = s c with _B[A(s c )] — <7[A(s c )]. Then we obtain by 
the help of (|B17|) that equation (|B12|) is not solvable at 



Next, we consider contributions to second-order vari- 
ational perturbation expansion ((H — -fftriai) 2 )triai,c ([53]) 
containing at least one factor cr a p. As described in the 
textbook Ref. HO for the case of the anharmonic oscil- 
lator, these terms can be best derived with the help of 
the square root trick. In our system this trick consists 
in substituting A(s) in A/ var of Eq. JTTJ) by (1 - k)A(s) 



denoted by A/ var (fc). The ({H-H tlial f 



' trial. c 



2 terms for 



k = containing at least one factor a a p are then given 
by 



-Htrial) ) trial, c, 2 



(B18) 



This leads to the contributions in A/ var 2 
(|B7p and the saddle point equation fld]) 
which are a factor (Ad/dA)g(A)/g(A) or 
(Ad/dA) 1+m g(A)/{Ad/dA)g(A) where m = 1,2 
smaller than the leading contributions. By using l|24[). 
(|27|1 we obtain only non-negligible contributions to / V ar,2 
(|B7j) or the saddle point equation (|B9jl . i.e. P t ^ 1, for 
A < 1/(Z ; (0) ) 4 / 3 in the fluid phase and A < l/(zi 1} ) 2 
in the solid phase. We point out that A max (|4T|) in the 
fluid phase and (|50|) in the solid phase in first-order 
variational perturbation theory is much larger than these 
A values. Note that we obtain also corrections Pi ^ 1 in 
the regime A ^> l/Z^ much larger than A max . 

Thus, we consider the regime A <C 1/(Z^) 4 ^ 3 in the 

fluid phase, and A <C l/(zj,) 2 in the solid phase. Here, 
we obtain prefactors Pi in (|B7p , (|B9|) which differ in gen- 
eral for both phases. We obtain 



_ 11 pT->oo _ 7 pT-^0 _ -i p 
— IT ) r \ — R ' r 2 ' L 



pT^O _ 11 pT^oo 
^3 — T ' 3 



iT— >oo 
2 



7 pT^O 

8 ' *4 



pT— too 
*4 



i . (B19) 
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